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ABSTRACT 

We report the detection of sub-mm emission from dust at 850 fim and of the 
^^CO J=4-3 hne in the two distant powerful radio galaxies 4C 60.07 (z=3.79) 
and 6C 1909+722 (z=3.53). In the case of 4C 60.07 the dust emission is also 
detected at 1.25 mm. The estimated molecular gas masses are large, of the 
order of - (0.5 - 1) x 10" Mq. The large FIR luminosities (Lpm ~ 10^^ 
suggest that we are witnessing two major starburst phenomena, while the 
observed large velocity widths (AVfwhm ^ 500 km s~^) are characteristic of 
mergers. In the case of 4C 60.07 the CO emission extends over ~ 30 kpc and 
spans a velocity range of ^ 1000 km s~^. It consists of two distinct features 
with FWHM of ^ 550 km s~^ and ~ 150 km and line centers separated by 
^ 700 km s~-^. The least massive of these components is probably very gas-rich 
with potentially ^ 60% of its dynamical mass in the form of molecular gas. 
The extraordinary morphology of the CO emission in this object suggests that 
it is not just a scaled-up version of a local Ultra Luminous Infrared Galaxy, 
and it may be a formative stage of the elliptical host of the residing radio-loud 
AGN. Finally we briefly explore the effects of the wide range of gas excitation 
conditions expected for starburst environments on the luminosity of high-J CO 
hues. We conclude that in unlensed objects, CO (J-l-1— > J), J-|-l>3 hues can be 
significantly weak with respect to CO J=l-0 and this can hinder their detection 
even in the presence of substantial molecular gas masses. 

Subject headings: galaxies: individual: 4C 60.07 — galaxies: individual: 6C 
1909-1-722 galaxies: interstellar matter — galaxies: formation — galaxies: active 
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1. Introduction 

The successful detection of CO emission in many local (z ^ 0.3) IRAS galaxies revealed 
large reservoirs of molecular gas mass that fuel nuclear starbursts and Active Nuclei 
(e.g. Tinney et al. 1988; Sanders, Scoville & Soifer 1991; Solomon 1997). In the most 
IR-luminous galaxies (LpiR ^ 10^^ M©) gas-rich mergers and interactions are thought to 
play a crucial role in causing the rapid accumulation of molecular gas in the center and thus 
initiating the onset of spectacular starbursts. Subsequent high resolution imaging of the CO 
emission in Ultra Luminous IR Galaxies (hereafter ULIRGs) (e.g. Scoville, Yun & Bryant 
1997; Downes & Solomon 1998) confirmed this picture by demonstrating the presence of 
molecular gas with high mass surface densities (S(H2) 5 x 10^ M© pc~^) and large 
velocity widths (~ 500 km s~^) confined within a few hundred parsec in the nuclear regions. 

Besides being responsible for some of the most spectacular starbursts in the local 
Universe the merging process is thought to play an important role in galaxy formation at 
high redshift, especially for spheroidal systems. The favorite formation model for these 
systems is that they grow from the hierarchical clustering of gas-rich "fragments" where the 
oldest stars (^ 10 Cyr) have already formed (e.g. Baron & White 1987; White 1996). The 
subsequent intense star formation rapidly consumes the gas mass and finally gives rise to a 
present-day giant-elliptical galaxy with its red colors, evolved stellar population and large 
mass, but devoid of substantial amounts of gas and dust. Powerful high-z radio galaxies 
(hereafter HzRGs) may be the progenitors of the massive present-day ellipticals hosting a 
radio-loud AGN. Since such galaxies seemed to have already settled into ellipticals (Best, 
Longair, & Rottgering 1997, 1998 and 1999) by z ~ 1, it is natural to assume that HzRGs 
at z > 2 with their frequently irregular morphologies (e.g. Pentericci et al. 1999) is where 
merger-induced large scale starbursts may have occured. 

The first detection of CO in IRAS 10214+4724 at z=2.286 (Brown & Vanden Bout 1991, 
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1992) and the subsequent detection of its mm/sub-mm continuum from dust (see Downes 
et al. 1992 and references therein) initiated ongoing efforts to detect CO and mm/sub-mm 
emission from the copious amounts of gas and dust expected in the high-z counterparts 
of ULIRGs and in galaxies undergoing their formative starbursts. The large negative 
K-corrections expected for the thermal dust spectrum and the high-J CO lines (e.g. Hughes 
1996; van der Wcrf & Israel 1996) of high-z galaxies as well as gravitational lensing help to 
make such objects detectable to the current mm/sub-mm instruments. Lensing amplifies 
the expected emission but usually complicates the gas and dust mass estimates. The 
galaxy IRAS 10214-1-4724 and many high-z CO/sub-mm luminous objects like the QSO 
H1413-I-117 (Barvainis et al. 1994), the recent sub-mm selected galaxies (Frayer et al 1998; 
Prayer et al. 1999) and the BAL quasar APM 08279+5255 (Downes et al. 1999) are lensed. 

Recently there have been detections of large amounts of dust in HzRGs (e.g. Hughes 
1996; Rottgering et al. 1998; Hughes & Dunlop 1998; Best et al. 1998), and a systematic 
survey is now been conducted to detect sub-mm emission from dust for all z > 3 radio 
galaxies (Rottgering et al. 1999) in order to understand the large range of FIR luminosities 
of these objects. The most notable examples are the extremely FIR-luminous radio galaxies 
8C 1435+635 (Ivison et al. 1998) and 4C 41.17 (Dunlop et al. 1994; Chini & Kriigel 
1994) whose properties do not appear to be influenced by gravitational lensing. In both 
cases the inferred gas masses are large enough (~ 10^^ M©) to suggest that the large FIR 
luminosities of these objects (~ 10^^ L©) are due to the formative starbursts. However, 
despite systematic efforts (Evans et al. 1996; van Ojik 1997) no CO emission has been 
detected from these two objects or any other powerful HzRGs. In this paper we present 
the detection of dust continuum and CO J=4-3 line emission in two powerful high-z radio 
galaxies and discuss the latter in the context of the earlier unsuccessful attempts to detect 
CO in this class of objects. Throughout this work we assume Ho = 75 km sec~^ Mpc~^ 
and qo = 0.5. 
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2. Observations and data reduction 

2.1. The SCUBA observations 

We used the Sub-mm Common User Bolometer Array (SCUBA) at the 15-m James 
Clerk Maxwell Telescope (JCMT)|^ in the photometry mode to observe 4C 60.07 on 1998 
April 15 and 6C 1909+722 on 1997 October 4 and 5 as part of an ongoing program to 
observe all HzRGs at z > 3 (Rottgering et al. 1999). SCUBA is a dual camera system 
cooled to ~ 0.1 K allowing sensitive observations with two arrays simultaneously. The 
short-wavelength array at 450 fim contains 91 pixels and the long-wavelength array at 
850 fim 37 pixels. The resolution is diffraction-limited with HPBW ~ 15" (850 /im), and 
HPBW ~ 8" (450 yum). For a full description of the instrument see Holland et al. (1998). 

We employed the recommended rapid beam switching at a frequency of 8 Hz and a 
beam throw of 60" in azimuth. The pointing and focus of the telescope were monitored 
frequently using CRL 618 and QSO 0836+710, and the typical rms pointing error was 
~ 3". Frequent skydips were used to deduce the atmospheric extinction. Typical opacities 
at 850 were r ~ 0.13 for our 1998 run and r ~ 0.20 for our 1997 run. Beam maps of 
Mars, CRL 618 and Uranus were used to derive the gain Cgso ~ 280 mJy/beam mV~^. The 
data were fiatfielded, corrected for extinction and sky noise was removed before they were 
co-added following standard procedures outlined in Stevens et al. (1997). The estimated 
rms of our measurements is consistent with what is expected from the total integration 
times and the NEFDs at the sky conditions of our runs. The systematic uncertainty in the 
flux density scale at 850 fim is of the order of ~ 10% (Holland, private communication). 

^The JCMT is operated by the Joint Astronomy Center in Hilo, Hawaii on behalf of 
the parent organizations PPARC in the United Kingdom, the National Research Council of 
Canada and the The Netherlands Organization for Scientific Research. 
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2.2. The IRAM interferometric observations 

We used the IRAM Plateau de Bure Interferometer (PdBI)[] between 1998 April 20 
and 1998 May 15 in the D configuration to observe CO J=4-3 (t'rcst = 461.040 GHz) in 
6C 1909+722 at z=3.534 (van Breugel private communication) in four tracks ranging from 
4-14 hours each. The same configuration was also used on 1998 November 8 and 29 to 
obtain two 8-hr tracks for 4C 60.07 at z=3.788 (Chambers et al. 1996). The correlator 
setup used for the 3 mm receivers involved 4 x 160 MHz modules covering a total bandwidth 
of 560 MHz tunned in single sideband. The band center was positioned at 96.290 GHz 
for 4C 60.07 and 101.770 GHz for 6C 1909+722. The 1 mm receivers were used in double 
sideband mode with the two remaining correlator modules as backends to simultaneously 
observe the continuum at 240 GHz. Typical SSB system temperatures were ~120 K 
(96.290 GHz) and ~150 K (101.770 GHz), and at 240 GHz they were 400 - 700 K (DSB). 
Bandpass calibration was obtained using 3C 454.3 and amplitude and phase calibration 
were obtained using IAP_0444+634 and NRAO 150. For the 3 mm receivers the residual 
phase noise was <^ 20° and for the amplitude <^ 5%. For the 1 mm receivers these figures 
are ^ 30° and ^ 20%. The flux density scale is accurate to within ~ 10%. 

After calibration and editing the data were processed with the standard NRAO AIPS 
software package. At 3 mm no continuum emission was detected in the line-free channels 
and no continuum subtraction was performed. The maps were produced using the task MX 
and in cases of low S/N, CLEAN was not applied. The rms noise of the maps is consistent 
with the thermal noise expected for the total observing time and average Tsys in the runs. 



^IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). 
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3. Results 

Both radio galaxies are detected in 850 fim and CO J=4-3 and are the brightest 
sub-mm objects found in the ongoing SCUBA survey of HzRGs at z > 3 so far, this being 
the main reason for selecting them for the follow-up sensitive CO observations. In 4C 60.07 
the CO (4-3) emission is clearly resolved (Figure 1). 

EDITOR: PLACE FIGURE |l| HERE. 

A remarkable characteristic of the CO emission in this galaxy is its large velocity 
range (^ 1000 km s~^) and its two distinct components with line centers separated by 
~ 700 km s~^. One component has AVfwhm ^ 550 km s^^ (extending beyond the observed 
band), while the other is narrower with AVfwhm ~ 150 km s^^. The feature with the 
narrow linewidth coincides with the position of the suspected radio core and thus the AGN, 
but a significant part of the broad linewidth component is clearly offset from it. This 
aspect of the CO J=4-3 emission is reminiscent of the CO J=5-4 emission detected in the 
radio-quiet quasar BR1202-0725 at z=4.69 (Omont et al. 1996) where a broad and a narrow 
linewidth component are also detected. In that object the narrow linewidth component is 
centered on the AGN position while the broad one is offset by ~ 4". 

The velocity-integrated map in Figure 2 shows extended CO emission with two peaks 
that are ~ 7" apart (~ 30 kpc at z = 3.791). Towards the same region 1.25 mm emission 
from dust is also detected, shown in Figure 3. Both the CO and 1.25 mm continuum maps 
are overlaid on a map of the non-thermal radio emission at 6 cm (Carilli et al. 1997). 

EDITOR: PLACE FIGURE | HERE. 



EDITOR: PLACE FIGURE | HERE. 
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The radio galaxy 4C 60.07 is a powerful radio galaxy (P4cm ~ 1-6 x 10^"^ W Hz~^) 
with a Fanaroff-Riley II (FR II) edge-brightened double lobe morphology (Fanaroff & Riley 
1974), hence there is the possibility of a significant contribution of the non-thermal emission 
to the observed 1.25 mm and 850 /im continuum. Extrapolating the non-thermal fiux 
density of Sgcm = 19 mJy of its brightest component with its associated spectral index of 
'^icm™ — (De Breuck, private communication) yields a negligible contribution 1%) 

in both 1.25 mm and 850 /im bands. Further confirmation of the thermal origin of the 
1.25 mm continuum is offered by the high resolution map in Figure 3 which shows that 
most of this emission is not associated with the observed non-thermal radio continuum 
and the peak 1.25 mm brightness is ~ 4" (17.6 kpc) offset from the radio core where the 
AGN probably resides (De Breuck private communication). Deep K'-band images obtained 
with Keck (van Breugel et al. 1998) reveal faint emission towards the weak eastern radio 
component but none from the regions with the brightest 1.25 mm continuum. Since at 
z ~ 3.8 the K'-band is rest-frame B-band, this is consistent with the brightest 1.25 mm 
emission marking the place of a massive gas/dust reservoir with large amounts of extinction. 

In order to obtain the highest S/N maps for the two distinct CO J=4-3 features we 
averaged all the appropriate channels and the two resulting maps are shown in Figure 4 
overlaid with the 1.25 mm continuum. From these maps it becomes obvious that the 
component with the largest dynamical and molecular gas mass is closely associated with 
the peak of the dust emission as expected. 

EDITOR: PLACE FIGURE | HERE. 

In the case of the radio galaxy 6C 1909+722 a map of the averaged CO J=4-3 emission 
and its spectrum are shown in Figure 5. There it can be seen that its linewidth is also 

^spectral index defined as S,^ oc z/" 
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rather large. Extrapolation of the non-thermal flux Sgcm = 59 mJy with the observed power 
law 020^ = — 1.3 (De Breuck private communication) yields a non-thermal contribution of 
~ 2% at 850 fim. All the observed parameters for the two HzRGs are tabulated in Table 1. 

EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE TABLE 1 HERE 

4. Discussion 

The two galaxies have extended radio emission and are unlikely candidates for 
gravitationally lensed objects. In the case of 4C 60.07, high-resolution radio images show 
a classical FR II source. Sensitive 6 cm (Carilli et al. 1997), K-band (Chambers et al. 
1996) and deep K -band images (van Breugel et al. 1998) do not reveal any obviously 
lensed features. For 6C 1909+722, HST images at 7000A do not reveal any lensed features 
either (Pentericci, private communication). Henceforth we assume that the two sources are 
unlensed and organize the discussion as follows: 

1. Estimate the molecular gas mass implied by the CO J=4-3 luminosity assuming 
global gas excitation conditions similar to the ones found in local starburst galaxies. 

2. Briefly explore the influence of various galactic environments on the molecular 
gas excitation conditions over large scales. Special emphasis is given on the effects of the 
various environments on the high- J CO lines and their detectability at high redshift. 

3. Use the mm/sub-mm data to flnd the dust masses and FIR luminosities and estimate 
the M(H2)/Mdust ratios, which we then compare to the values in the local Universe. 

4. Discuss the evolutionary status of these objects in terms of their gas mass relative to 
their dynamic mass, and their star formation rates and efficiencies. We focus particularly 



on 4C 60.07 where the two distinct molecular gas components strongly suggest a large scale 
starburst event unhke the usually nuclear starbursts observed in local ULIRGs. 



4.1. Molecular gas content 

The estimate of molecular gas mass from the CO J=l-0 luminosity involves the 
so-called standard conversion factor Xco (e.g. Young & Scoville 1982; Bloemen 1985; 
Dickman et al. 1986; Young & Scoville 1991), which is "calibrated" using molecular clouds 
in the Milky Way. Prom the CO(4-3) luminosity M(H2) is then given by 



M(H2 



X 



CO 



143 ^ 2 k l^? 



43 



1 + z 



Av 



K dv, 



where Dl = 2 c H~-'^(l + z — \/l + z) is the luminosity distance for Qo = 0.5, 2/43 is the 
rest-frame frequency of the CO J=4-3 transition, r43 is the (4-3)/(l-0) line ratio of the 
area/velocity-integrated brightness temperatures and '^v^hB observed flux density. 

Substituting astrophysical units yields 



a/Xm\ (1 + z - VI + z) 



/av S^obsdv 
Jy km s~^ 



(2) 



The dependence of Xco on the ambient conditions of the H2 gas has been extensively 
explored (e.g. Bryant & Scoville 1996; Sakamoto 1996). An important recent result is that 
in the intense starburst environments of ULIRGs the warm, diffuse gas phase dominating 
the ^^CO emission does not consist of virialized clouds, hence leading to an overestimate 
of M(H2) when the standard Xco ^ 5 (km s~^ pc^)~^ is used (Solomon et al. 1997; 
Downes & Solomon 1998). A more appropriate value for such ISM environments and hence 
their high-z counterparts is Xco ~ 1 (km s""*^ pc^ (Downes & Solomon 1998). 
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A plausible value for 143 can be found assuming an ISM environment similar to that 
of a local "average" star burst. The average hne ratios measured towards such systems 
are 121 ~ 0.9 (e.g. Braine & Combes 1992; Aalto et al. 1995), 1^2 ~ 0.64 (Devereux 
et al. 1994), while the ^^CO/^^CO J=l-0, 2-1 ratios are Rio ^ R21 ^ 13 (Aalto et 
al. 1995). A Large Velocity Gradient (LVG) code with the aforementioned ratios as 
inputs was employed to model the physical conditions of the gas (e.g. Richardson 1985). 
The conditions corresponding to the best fit are Tkm = 50 K, n(H2) ~ 10^ cm"*^ and 
[C0/H2]/dV/dr = 3 x 10~^ (km s~^ pc"^)"^, similar to the ones deduced for the high-z 
starburst IRAS 10214+4724 (Solomon, Downes, Radford 1992). For these conditions it is 
r43 = 0.40, which for Tcmb = (1 + z) x 2.75 K ~ 13 K, is slightly enhanced to r43 = 0.45. 
This ratio and Xco = 1 (km s'^ pc^)-^ yield M(H2) ^8x lO^o Mq for 4C 60.07 
and M(H2) ^ 4.5 x 10^° M© for 6C 1909+722. These are most hkely lower hmits since 
even in ULIRGs the ratio r43 can be significantly smaller but not much larger. Indeed the 
assumption of an "average" starburst excitation environment, while convenient, hides the 
fact that there is a wide range of physical conditions for the gas in IR-luminous galaxies, 
and in much of that range C0(J+1— >J), J + 1 > 3 can be significantly weaker than CO (1-0). 

As we shall see, the potential faintness of the high-J CO transitions may partly explain 
why attempts to detect high-z CO have been less successful than the ones trying to detect 
mm/sub-mm continuum from dust. Especially for HzRGs two large systematic searches 
(Evans et al. 1996; van Ojik et al. 1997) gave null results. On the other hand observations 
of the sub-mm continuum fare better with 7 such objects detected at 850 //m (see van der 
Werf 1999 and references therein). A similar situation applies also in other types of high-z 
objects. While the uncertainties of the gas-to-dust ratio and a rising Lfir/Lco ratio with 
FIR luminosity (van der Werf 1999) can still adequately account for this, the wide range of 
gas excitation observed in IR-luminous galaxies cannot but have an effect on the luminosity 
of CO lines especially at high J levels. 
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4.2. The high-J CO transitions in ULIRGs at high z 

In the attempts to detect CO in high-z objects various transitions are observed 
depending on the particular redshift and the receiver used. For z ;^ 2 mostly C0(J+1^J) 
with J + 1 > 3 is observed and in the systematic searches in HzRGs transitions with 
J+l=4-9 were routinely observed (Evans et al. 1996; van Ojik 1997). It has been argued 
that observing higher J transitions can offset the dimming due to the distance in high-z 
objects in a manner similar to the negative K-corrcctions of the thermal spectrum from 
dust (van der Werf & Israel 1996). In this picture the warm and dense gas in an "average" 
starburst environment thermalizes the CO transitions up to J-|-l=6. However the luminosity 
of the high-J CO transitions may be much smaller for two basic reasons, namely 

1) The presence of a warm (Tkin = 50 - 100 K) but diffuse (n(H2) ~ 10^ - 10^ cm-^) 
and subthermally excited (for J -|- 1 > 2) gas phase that dominates the ■'^^CO emission in 
ULIRGs (Aalto et al. 1995; Downes & Solomon 1998). The most conspicuous such galaxy is 
Arp 220, frequently used as the standard ULIRG for comparison with high-z FIR-luminous 
sources. In this source a low r2i = 0.53 and high Rio > 20, R21 = 18 ratios are observed 
(Aalto et al. 1995; Papadopoulos & Seaquist 1998), in contrast to the "average" starburst 
ratios mentioned previously. Another important characteristic of this gas phase is the 
moderate optical depths (r ~ 1 — 2) of the ^^CO (1-0) transition (Aalto et al. 1995). 

2) A large reservoir of cold and/or subthermally excited gas extending beyond the 
warm star bursting nuclear region of an ULIRG. Such excitation gradients are observed 
when large beams that sample the extended emission are used (Papadopoulos & Seaquist 
1998). The associated cool dust in the ULIRG VV 114 was imaged recently with SCUBA 
(Prayer et al. 1999) and extremely cold gas (Tkin = 7 — 10 K) is inferred over large scales 
for starbursts hke IC 5135 and NGC 7469 (Papadopoulos & Seaquist 1998). This gas phase, 
if present, can easily dominate the global CO excitation especially in high-z systems where 
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a beam of ;^ 5" at z ;^ 2 corresponds to L ;^ 30 kpc. 

Here we must stress that the effectiveness of a cold gas component in suppressing the 
observed global gas excitation is not altered by the higher CMB temperature at high z. 
Indeed, while the higher CMB temperature enhances the populations of the high J levels 
of CO, it also corresponds to a higher background against which the respective lines must 
be detected. Moreover the effective temperature of cold gas at high z is not simply the 
sum of the temperature of this gas phase at z—0 and (1 + z) x 2.75 K as one might naively 
assume. This can be demonstrated with some simple arguments. Obviously the excitation 
temperature Texc of any coUisionally excited line at any redshift is bounded as 



where Tcmb — 2.75 K is the present epoch CMB temperature and Tkin is the gas 
kinetic temperature. If we further assume that on large scales in the ISM of a galaxy 

thermodynamic equilibrium exists between gas and dust, i.e. Tjust = ^km, then for a dust 
emissivity of a = 2 the energy balance within a typical giant molecular cloud is described by 



where Uisrf is the energy density of the interstellar radiation field (O, B, A stars) that 
heats the grains directly and Umech is the "mechanical" heating energy density deposited 
to the molecular clouds (e.g. supernovae shocks, turbulent cascade, cloud-cloud collisions) 
that first heats the gas and then the dust. These are the two major heating mechanisms 
of an average molecular cloud and since, a) more cooling processes are at play than only 
dust radiation (e.g. C"*", CI, CO hues), b) usually Tgas < T^ust (except in the surfaces of 
UV-illuminated clouds, HoUenbach & Tielens 1997), it follows that the dust temperature 



(1 + Z) X Tcmb < Texc < Tkin, 



(3) 



Uisrf + U^ech = cooling rate oc Tdust(z)^ - 2.75^ x (1 + zf , 



(4) 
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yielded by the last equation will be an upper limit to the gas temperature. 

The aforementioned physical processes responsible for the heating of the ISM do 
not depend on the particular value of the CMB temperature at a given z. Thus the last 
equation, being redshift-invariant, allows us to find the high-z equivalent temperature of any 
given ISM component from its temperature in the local Universe (see also Combes, Maoli 
& Omont 1999). It is easy to see that at high z, the low-temperature component of the 
the local ISM has T^ust (1 + ^i) x Tcmh which also means (Equation 3) that lines become 
thermalized at the expense of becoming invisible against the enhanced CMB radiation field. 

Choosing a nominal redshift of z = 3.5 as representative of the redshifts of the 
two HzRGs we find that for a warm ISM component of Tdust(z = 0) = 50 K the 
Tdust(z — 3.5) is essentially identical. However for a cold component of T^ust = 10 K it 
is Tdust(z = 3.5) = 12.892 K, i.e. just ~ 0.5 K above the (1 + z) x 2.75 K = 12.375 K 
temperature of the CMB at that redshift. In Table 2 we display the expected hne ratios at 
z = 3.5 for three sets of conditions, namely a) an "average" starburst, b) diffuse and warm 
gas of moderate ^^CO J=l-0 optical depth, and c) cold gas, representative of the phase 
that may exist in ULIRGs beyond the central starburst region. 

EDITOR: PLACE TABLE 2 HERE 

This table shows that a wide range is expected for CO(J-l-l^J), J>2 hne luminosities 
relative to CO J=l-0. Several published r32 = (3—2) /(I — 0) line ratios (Lisenfeld et al. 1996) 
that sample the global CO emission in ULIRGs, and a recent survey of this ratio in many 
nearby galaxies (Mauersberger et al. 1999) reveal a large range of values, namely 
r32 ~ 0.1 — 1, thus a similar or larger range is to be expected for the higher J transitions. 
Moreover since the flux density ratio S(J + 1 — J)/S(l — 0) = (J + 1)^ rj+i j determines 
whether a high- J transition is easier to detect than J=l-0 (assuming equal sensitivities), it 
can be seen (Table 2) that CO J=3-2 is the highest transition for which this ratio is ;^ 1 
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for all the expected conditions. Hence many non-detections of C0(J-|-1— >J), J>2 in high-z 
systems can be due to the globally sub-thermal excitation of these lines rather than a true 
molecular gas mass deficiency. At the same time this makes the upper limits for M(H2) from 
such non-detections much less stringent. For example if the CO J=5-4 hne was observed, 
and for Xco = 1 (km s"^ pc'^)'^, the upper limit in M(H2) can be ~ 1/0.025 x 1/5 = 8 
times larger that what is usually reported for this line (e.g. van Ojik 1998). 

Finally for gravitationally lensed objects with underlying steep gas excitation gradients, 
differential amplification can render line ratios useless in deducing the average gas excitation 
unless an accurate description of the lensing potential is available. Such gradients are 
expected in nuclear starbursts where gas warm and dense enough to excite high-J CO 
lines is more centrally confined in the immediate area of the starburst. Recent mapping 
of CO J=4-3 in M 51 and NGC 6946 (Nieten et al. 1999) showed that even in more 
quiescent galaxies the highly excited gas is strongly concentrated in the center. Under 
such circumstances differential amplification can enlarge the effective size of the high-J 
CO emitting region with respect to the low-J one and thus alter the observed global line 
ratios towards the ones expected for warm and dense gas. The frustrating aspect of this 
effect is that it can become progressively severe for lines that are widely separated in J, 
i.e. exactly the ones whose ratios are most sensitive to the gas excitation conditions under 
normal circumstances. 

Clearly more observations of high-J CO fines in local ULIRGs are needed in order to 
reveal the range of their global luminosities and relative brightness distributions and thus 
allow a better understanding of similar systems at high-z. 
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4.3. Dust mciss, the geis-to-dust ratio and LpiR 

There are many uncertainties associated with the estimate of the dust mass from 
mm/sub-mm measurements (e.g. Gordon 1995), with the uncertainty of the dust 
temperature being one of the most important ones (Hughes 1996). In the case of 4C 60.07 
the detection of the dust continuum both at 1.25 mm and 850 /im allows us to place some 
broad constraints on the dust temperature as long as the emissivity law is a > 1. Indeed, 
for an optically thin isothermal reservoir of dust mass it is 

R(a, Td) = = 1.47" ^ , (5) 

^ ^ Si.25mm (e^VTd - 1)-^ - 0.0162 ' 

where is the dust temperature and a = 1 — 2 is the exponent of the emissivity law. 

For 4C 60.07 we find R(a, Td) = 2.44 ± 0.73. We adopt an emissivity law of a = 2, 
which was deduced in other high-z objects, e.g. IRAS 10214+4724 (Downes et al. 1992) 
and 8C 1435+635 (Ivison et al. 1998) over a similar rest-frame spectral range. For this 
emissivity law the ilcr range of the observed R{a, Td) yields a range of Td ~ 20 — 50 K. 
Given the starburst nature of these two objects and their luminous CO (4-3) line (the J=4 
level is ~ 55 K above the ground state) we assume a dust temperature of Td = 50 K. Then 
the dust mass can be estimated from 



Mdust = ^\ ^'^7 , [B(i^em, Td) - B(i/em, Te^b(z))]~' , (6) 

(1 + Z) kd(z/em) 

where i^em = (1 + ^)^obs is the emitted frequency, B(i/, T) is the Planck function, Tcinb(z) 
is the CMB temperature at redshift z, and kd(i^em) = 0.04(i/em/250GHz)^ m^ kgr~^ is the 
adopted dust emissivity (e.g. Kriigel, Steppe, & Chini 1990). For the assumed Td the 
CMB term can be omitted with ^ 3% error. Here it is worth noting that for cold dust 
this term can be significant. For example Tcmblz = 4) = 13.75 K, hence a dust component 
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with Td = 15 K at that redshift is a very cold component, just ~ 1 K above Tcinb(z = 4). 
An observed wavelength of 1.25 mm corresponds to hz/em/k ~ 57.4 K at z = 4 and for this 
wavelength the CMB term is ~ 70% of the dust term. 

Assuming similar parameters for the dust emission in 6C 1909+722 the 850/xm flux 
density yields M^ust = 1.5 x 10^ Mq for both objects. This gives warm gas-to-dust ratios 
of M(H2)/Mdust ~ 530 (4C 60.07) and M(H2)/Mdust ~ 300 (6C 1909+722), i.e. within the 
range of the values found for local IRAS galaxies (e.g. Young et al. 1986; Stark et al. 1986; 
Young et al. 1989) and ULlRGs (Sanders et al. 1991). This suggests that in a Universe 
at ~ 10% of its current age these two HzRGs already have heavy-element abundance 
comparable to the galaxies in the contemporary Universe. It is also worth noting that, while 
many uncertain factors enter into the estimate of the gas-to-dust ratio, the values chosen 
for them here are the ones considered plausible from the study of the ISM in local ULIRGs. 

The FIR luminosities of the two HzRGs were estimated by assuming the underlying 
spectrum of an optically thin, isothermal reservoir of dust mass, hence 

LfIR= / Lr/emdi^em = 47rMdust / kd(i/em)B(l/em, Td)dl/em, (7) 

Jo Jo 
after using Equation 6 to substitute Mdust this finally yields 

Lfir = 47rA(a) x-("+^) (e^ - 1) S,„,^i/obs, (8) 

where x — hi/em/kTd and X{a) a numerical constant that depends on the emissivity law. In 
astrophysical units this becomes 

LpiK = 2 X 10^ (l + z - VTT^fxia) x-("+^) (e^ - 1) (^) L^. (9) 
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For the assumed T^ust = 50 K the last expression will yield a lower hmit to the true LpiR 
since in some ULIRGs the optical depth may become significant even at FIR frequencies 
(e.g. Solomon et al. 1997), and in a starburst environment dust can be warmer still. We 
find LpiR 1.5 X 10^^ L0 for both HzRGs, a luminosity comparable to 8C 1435+635 and 
4C 41.17 (Ivison et al. 1998 and references therein) the other two high-z radio galaxies 
detected in the sub-mm (rest-frame FIR) spectral range whose properties do not appear to 
be infiuenced by gravitational lensing. 

4.4. The evolutionary status of 4C 60.07 and 6C 1909+722 

The usual conclusion drawn when such large FIR luminosities are deduced for a galaxy 
is that we are witnessing it during a starburst phase. However many of the IR-luminous 
galaxies in the local Universe and certainly the two particular HzRGs also harbor an AGN, 
thus part of the FIR luminosity maybe due to AGN-heated dust. Nevertheless, recent 
studies of ULIRGs (Genzel et al. 1998; Downes & Solomon 1998) reveal that most of them 
are powered by recently formed massive stars. Also, even when an AGN does make a 
significant contribution to Lfir, it is usually ^ 30% (Genzel et al. 1998; Downes & Solomon 
1998). Therefore the quahtative arguments behind converting FIR luminosities to star 
formation rates are not hkely to change even if an AGN is present. For 4C 60.07 the case 
for the starburst origin of its FIR luminosity is made stronger by the fact that most of its 
1.25 mm dust emission is not emanating in the vicinity of the AGN. 

The FIR luminosity provides a measure of the current star formation rate (SFR) 
according to 

SFR = 10-^° *LpiR Mo yr-S (10) 
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where ^ ~ 1 — 6, depending on the assumed IMF (Telesco 1988 and references therein). 

Adopting the most conservative value ^ ~ 1, the FIR luminosity of the two HzRGs 
imphes SFR ~ 1500 yr~^. Such a high rate of star formation can, if sustained, produce 
10^^ — 10^^ Mq of stars in ~ 0.06 — 0.6 Gyr and it is comparable to the SFRs found for 
8C 1435+635 and 4C 41.17 (Ivison et al. 1998). The efficiency with which such a burst 
converts molecular gas into stars is given by SFR/M(H2) or equivalently by Lfir/M(H2) 
(assuming the same ^ for all the galaxies). On many occasions the quantity Lfir/Lco 
is used instead, but since the Xco conversion factor is ~ 4 — 5 times smaller in ULIRG 
systems the latter ratio can be misleading when a single value for Xco is used for galaxies 
spanning a FIR range from moderately IR-luminous galaxies (Lpm ~ lO^'^ L©) to ULIRGs 

(Lfir;^io12 Lo). 

The star formation efficiencies estimated are ~ 190 M^^ (4C 60.07) and 
~ 330 L0 M"^ (6C 1909+722), comparable to the ones found for ULIRGs (Solomon et 
al. 1997) once M(H2) has been estimated using the same Xco factor. The implied star 
formation rates and efficiencies of these galaxies point towards spectacular starbursts 
occurring at high redshift. In addition the large linewidths observed, which in the case of 
4C 60.07 exceed ~ 1000 km s~^, are routinely observed towards ULIRGs (e.g. Solomon et 
al. 1997) and are a kinematic signature of the ongoing mergers/interactions that trigger 
their enormous starbursts (Sanders, Scoville, & Soifer 1991). 

The tantalizing question raised when high-z starbursts are found is to what extent 
the observed star formation episode is forming the bulk of their eventual stellar mass. 
The main difficulty in answering this important question stems from the uncertainties 
involved in deducing the total gas mass from CO measurements. If we assume an 
atomic-to-molecular gas ratio of M(HI)/M(H2) ~ 2 found for IRAS galaxies (Andreani, 
Casoh, & Gerin 1995) and a large sample of spirals (Casoli et al. 1998), we obtain a total 
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gas content of Mga^ = 1.0 x 10^i(XcoA43) (4C 60.07) and Mg^^ = 6.0 x 10^°(XcoA43) 
(6C 1909+722). For the "average" starburst value = 0.45 this gives gas masses of the 
order of ~ (1.3 — 2.2) x 10^^ Mq. These are an order of magnitude higher than typical gas 
masses in local ULIRGs (e.g. Downes and Solomon 1998) and constitute ~ (15 — 20)% of 
the total stellar mass of a typical elliptical associated with the 3CR radio galaxies at z ~ 1 
(Best, Longair & Rottgering 1998). Thus it is clear that these enormous starbursts still 
have a vast reservoir of gas to be eventually turned into stars. 

The best evidence yet that at least in the case of 4C 60.07 we are witnessing an 
extraordinary starburst and not a scaled-up version of an ULIRG comes from the fact that 
in this galaxy the CO J=4-3 emission from gas is distributed over scales of ~ 30 kpc. This 
is in contrast with the local ULIRGs where most of the molecular gas and starburst activity 
is confined within the central ^1 — 5 kpc. Since starburst activity is usually co-extensive 
with the molecular gas reservoir, and this will be particularly true for the high excitation 
CO J=4-3 hne, it seems that the starburst in 4C 60.07 occurs on galaxy-wide scales. 

Equally intriguing is the fact that the CO emission consists of two distinct components 
widely separated in velocity (Figures 1, 4). Their velocity- integrated CO flux densities of 
^ 1.65 ± 0.35 Jy km s^^ (wide-linewidth component) and 0.85 ± 0.20 Jy km s"^ (narrow- 
hnewidth component) yield masses of M(H2) ;^ 5 x 10^° M© and M(H2) ^ 2.6 x 10^° Mq 
respectively (Equation 2, r43 = 0.45). We derived size estimates for these two components 
by fitting both the image and the visibility plane of the CO emission shown in Figure 4 with 
an underlying gaussian brightness distribution. The narrow-linewidth component appears 
unresolved with an upper limit of ~ 4" while the wide-linewidth one is marginally resolved 
with the largest size being ~ 5". In the absence of adequate spatial/kinematic information 
allowing a distinction between the possible geometrical arrangements of the CO-cmitting 
gas we assume that the largest estimated size L corresponds to the diameter of a disk, hence 
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its mass is given by 

- ^il^ - "'^ X ^"^ (ife)^ (^) 0"' (U) 

where i is the inchnation of the disk and aj ~ 1 (Bryant & Scoville 1996). 

For the wide-hnewidth component the largest estimated size corresponds to L ~ 22 kpc, 
which gives Mdyn ~ 7.7 x 10-'^^(sin^ i)~^ Mq, comparable to the mass of a present-day 
giant elliptical. The ratio of the inferred molecular to dynamic gas mass then is 
M(H2)/M(jyn ~ 0.06 sin^ i. For the narrow-velocity component the upper limit corresponds 
to L ~ 17.5 kpc yielding a dynamic mass of Mdyn ^ 4.6 x 10^°(sin^ i)"-*^ Mq. This is 
^ 5% of the mass of a typical giant elliptical and a comparison with its molecular gas 
mass gives M(H2)/Mdyn ^ 0.60 sin^ i. Thus, geometrical factors aside, this component is 
significantly richer in molecular gas than the more massive one and it maybe the one where 
star formation has yet to form the bulk of its eventual stellar mass. 

It is easy to contrive a combination of different excitation properties, Xco values, 
and inclinations so that M(H2)/Mdyn in one or both molecular gas components is altered 
significantly. Still the fact that remains unaltered is that in 4C 60.07, unlike in a typical 
ULIRG, the intense star formation occurs over large scales in two spatially and kinematically 
distinct molecular gas reservoirs. This brings in mind the scenario for the formation of a 
giant elliptical at high redshift where several star-forming clumps merge to eventually form 
the galaxy. In this picture a gas-rich low mass clump like the one seen in 4C 60.07 is still in 
the process of merging and vigorous star formation while the higher mass object has already 
formed most of its stars. Of course the possibility that the wide-linewidth component itself 
consists of several virialized gas-rich clumps that do not necessarily constitute a bound 
system cannot be discarded in the light of the present data. 
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For 6C 1909+722 the CO emission is unresolved with a size ^ 4" (L ^ 18 kpc), which 
for the observed hnewidth (Table 1) corresponds to Mdyn ^ 5.85 x 10^^(sin^ i)"^ Mq and 
M(H2)/Mdyn ^ 0.08 sin^ i. Thus, while there is little doubt that this galaxy is a starburst, it 
can still be a high-z counterpart of an ULIRG where most of its eventual stellar population 
has already formed and the intense star formation is confined in its central region. 

Finally, the case of 4C 60.07 makes it clear that when searching for high-z CO lines a 
large velocity coverage is necessary not only because of their often uncertain redshift, but 
in order to be able to detect various components that may be far apart in velocity. Higher 
resolution observations are necessary in order to reveal more structure in the molecular gas 
reservoir of this galaxy and allow better constraints of its dynamical mass. However the 
faintness of the CO emission observed here may hinder such efforts and a systematic study 
of such objects may have to await the advent of the next generation mm interferometers 
with significantly larger collecting areas. 

5. Conclusions 

In this paper the detection of mm/sub- mm emission from dust and the first detection 
of CO J=4-3 in two powerful high-z radio galaxies has been presented. The analysis of the 
data leads us to the following conclusions: 

1. Using the most conservative CO/H2 conversion factor (1/5 of the galactic value) the 
CO J=4-3 emission implies molecular gas masses of the order of M(H2) ~ (0.5 — 1) x 10^^ 

in the two high-z powerful radio galaxies. Their sub-mm continuum corresponds to large 
FIR luminosities (~ 10^^ Lq) which imply that we are witnessing intense starburst 
phenomena at z ~ 3.5 — 3.8 converting the aforementioned gas mass into stars. 

2. The wide range of the molecular gas excitation expected in Ultra Luminous IR 
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galaxies over large scales is briefly explored and we conclude that observing objects with 
similar ISM conditions at high z using high-J CO transitions may in some cases hinder their 
detection. This can partly explain explain why, despite systematic efforts, CO detections 
of unlensed objects are fewer than detections of their mm/sub-mm continuum from dust. 
Therefore it is still too early to draw any general conclusions about the abundance of 
molecular gas in high-z systems from the data currently available in the literature. 

3. The estimated molecular gas-to-dust ratios found in these two objects are in the 
range found for the local IRAS galaxies, revealing that in both of them most of the heavy 
elements have already been produced at present-day abundances. 

4. In 4C 60.07 the CO emission and presumably the starburst activity implied by its 
large FIR luminosity is distributed over scales of ~ 30 kpc and consists of two distinct 
components spanning a total velocity range of ^ 1000 km s~^. Thus this galaxy does not 
seem to be just a scaled-up high-z version of a local Ultra Luminous IR galaxy where 
most of the starburst activity and the accompanying bright CO/dust emission are confined 
to a more compact region in the center. A plausible scenario is that we are witnessing 
the ongoing formation event of a giant elliptical galaxy, the future host of the residing 
radio-loud AGN. 
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Fig. 1. — Naturally weighted channel maps of CO J=4-3 for 4C 60.07 at a spatial resolution 
of 8.9" X 5.5" and frequency resolution of 35 MHz (~ 110 km s"-^). The contours are 
(—3, —2, 2, 3, 4, 5, 6) x cr, with cr = 0.55 mJy beam~^. The FWHM of the restoring beam is 
shown at the bottom left of the first channel. The cross marks the position of the radio core, 
RA (J2000): 05 12 55.147, Dec (J2000): +60 30 51.0 (De Breuck private communication) 

Fig. 2. — Uniform weighted, velocity-integrated map of CO J=4-3 in 4C 60.07 (contours) 
overlaid on a 6 cm map (greyscale) at a common resolution of 7.1" x 5.2" The contours are 
(2,3,4,5) X cr, with cr = 0.30 Jy beam~^ km s~^ and the greyscale is: 1.5-15 mJy beam~^. 
The FWHM of the beam is shown at the bottom left. Note: CLEAN was not applied in the 
CO J=4-3 uniform- weighted channel maps because of the low S/N in individual channels. 
The cross marks the position of the radio core (see Figure 1). 

Fig. 3. — Naturally weighted map of the 1.25 mm continuum emission of 4C 60.07 (contours) 
overlaid to 6 cm continuum emission (greyscale) at a common resolution of 3.7" x 1.9". The 
contours are (—3,-2,2,3,4,5) x cr, with cr = 0.6 mJy beam~^ and the greyscale range is 
0.6-12 mJy beam~^. The FWHM of the restoring beam is shown at the bottom left. The 
cross marks the position of the radio core (see Figure 1). 

Fig. 4. — Naturally weighted CO J=4-3 maps of the two kinematic components of 4C 60.07 
(contours) overlaid to its 1.25 mm continuum (greyscale), at a common resolution of 
8.9" X 5.5". The velocity offset is denoted at the upper right, and the FWHM of the line at 
the bottom right. The contours are (—3, —2, 2, 3, 4, 5, 6, 7) x cr, with cr = 0.3 mJy beam~^ 
(top) and a = 0.4 mJy beam~^ (bottom), and the greyscale range is 1.3-5.85 mJy beam~^. 
The FWHM of the CLEAN beam is shown at the bottom left. The cross marks the position 
of the radio core (see Figure 1). 
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Fig. 5.— The CO J=4-3 spectrum of 6C 1909+722 (top) with Ai/^han = 50 MHz 
(~ 145 km s~^) and naturally weighted maps of the CO J=4-3 emission centered at 
V = 101.725 GHz (middle) and u = 101.975 GHz (bottom). The averaged frequency 
interval in both maps is Au — 200 MHz (~ 590 km s~^) and the contours are 
(—3, —2, 2, 3, 4, 5, 6, 7, 8, 9) x cr, with a — 0.3 mJy beam~^. The cross marks the map center 
(see Table 1) which is within 0.5" of the position of the radio core (De Breuck, private 
communication). The beam size is 7.65" x 7.11" and its FWHM is shown in the bottom 
map. 



-33- 
Table 1 

Observational parameters 
Parameter 4C 60.07 6C 1909+722 



RA (J2000)^ 
Dec (J2000)^ 

2! CO 

Ico (Jy km s"^) 
AVfwhm (km s"^) 

Sl.25mm (mjy) 

Ssmm (mJy) 
SssoMm (mJy) 



05^^ 12™ 54^80 
+60° 30' 51.7" 
3.791 
2.50 + 0.43 

;^iooo 

4.5 + 1.2 

< 0.5 {2a) 
11.0 + 1.5 



19^^ 08^^ 23^70 
+72° 20' 11.8" 
3.532 
1.62 + 0.30 
530 + 70 
< 2 (2(7) 
< 0.6 (2a) 
13.5 + 2.8 



Coordinates of the image center. 



Table 2 



ISM environments and line ratios at z = 3.5 



Line ratios 


"Average" starburst 


Diffuse warm gas 


Cold gas 


J=2-l 


0.95-1.12 


0.64-1.06 


0.40-0.61 


J=3-2 


0.73-1.07 


0.30-0.92 


0.14-0.29 


J=4-3 


0.45-1.04 


0.10-0.74 


0.035-0.11 


J=5-4 


0.21-0.99 


0.025-0.52 


0.005-0.02 


J=6-5 


0.065-0.93 


0.006-0.31 


< 10-3-0.002 


J=7-6 


0.015-0.84 


0.0015-0.15 


< 10-3 


J=8-7 


0.002-0.69 


< 10-3-0.05 


< 10-3 


J=9-8 


< 10-3-0.50 


< 10-3-0.013 


< 10-3 



Note. — All tlie line ratios are normalized by the ^^CO J=l-0 brightness, and the 
assumed CMB temperature is (1 + 2;) x 2.75 K = 12.375 K. 

"Average" starburst: Tkm = 50 - 100 K, n(H2) = 103 - 10^ cm-3, 
[C0/H2]/(dV/dr) = 3 X 10-^ (km s-^ pc-^)-^ 
Diffuse warm gas: Tki„ = 50 - 100 K, n(H2) = 10^ - 103 cm-3, 
r ~ 1 - 2 for J=l-0. 

Cold gas: Tkin - 13 K (see text), n(H2) = 10^ - 103 cm-3, 
[C0/H2]/(dV/dr) = 3 X 10-^ (km s-^ pc-^)-^ 
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